Automated Pharmacogenomic Reports for Clinical Genome Sequencing.

Clinical laboratories offering genome sequencing have the opportunity to return pharmacogenomic findings to patients, providing the added benefit of preemptive testing that could help inform medication selection or dosing throughout the lifespan. Implementation of pharmacogenomic reporting must address several challenges, including inherent limitations in short-read genome sequencing methods, gene and variant selection, standardization of genotype and phenotype nomenclature, and choice of guidelines and drugs to report. An automated pipeline, lmPGX, was developed as an end-to-end solution that produces two versions of a pharmacogenomic report, presenting either Clinical Pharmacogenetics Implementation Consortium or US Food and Drug Administration guidelines for 12 genes. The pipeline was validated for performance using reference samples and pharmacogenetic data from the Genetic Testing Reference Materials Coordination Program. To determine performance and limitations, lmPGX was compared with three additional publicly available pharmacogenomic pipelines. The lmPGX pipeline offers clinical laboratories an opportunity for seamless integration of pharmacogenomic results with genome reporting.

Pharmacogenomics: novel loci identification via integrating gene differential analysis and eQTL analysis.

Nearly one-half of asthmatic patients do not respond to the most commonly prescribed controller therapy, inhaled corticosteroids (ICS). We conducted an expression quantitative trait loci (eQTL) analysis using >300 expression microarrays (from 117 lymphoblastoid cell lines) in corticosteroid (dexamethasone) treated and untreated cells derived from asthmatic subjects in the Childhood Asthma Management Program (CAMP) clinical trial. We then tested the associations of eQTL with longitudinal change in airway responsiveness to methacholine (LnPC20) on ICS. We identified 2484 cis-eQTL affecting 767 genes following dexamethasone treatment. A significant over-representation of lnPC20-associated cis-eQTL [190 single-nucleotide polymorphisms (SNPs)] among differentially expressed genes (odds ratio = 1.76, 95% confidence interval: 1.35-2.29) was noted in CAMP Caucasians. Forty-six of these 190 clinical associations were replicated in CAMP African Americans, including seven SNPs near six genes meeting criteria for genome-wide significance (P < 2 × 10(-7)). Notably, the majority of genome-wide findings would not have been uncovered via analysis of untreated samples. These results indicate that identifying eQTL after relevant environmental perturbation enables identification of true pharmacogenetic variants.

Genome-wide association analysis in asthma subjects identifies SPATS2L as a novel bronchodilator response gene.

Bronchodilator response (BDR) is an important asthma phenotype that measures reversibility of airway obstruction by comparing lung function (i.e. FEV(1)) before and after the administration of a short-acting ß(2)-agonist, the most common rescue medications used for the treatment of asthma. BDR also serves as a test of ß(2)-agonist efficacy. BDR is a complex trait that is partly under genetic control. A genome-wide association study (GWAS) of BDR, quantified as percent change in baseline FEV(1) after administration of a ß(2)-agonist, was performed with 1,644 non-Hispanic white asthmatic subjects from six drug clinical trials: CAMP, LOCCS, LODO, a medication trial conducted by Sepracor, CARE, and ACRN. Data for 469,884 single-nucleotide polymorphisms (SNPs) were used to measure the association of SNPs with BDR using a linear regression model, while adjusting for age, sex, and height. Replication of primary P-values was attempted in 501 white subjects from SARP and 550 white subjects from DAG. Experimental evidence supporting the top gene was obtained via siRNA knockdown and Western blotting analyses. The lowest overall combined P-value was 9.7E-07 for SNP rs295137, near the SPATS2L gene. Among subjects in the primary analysis, those with rs295137 TT genotype had a median BDR of 16.0 (IQR = [6.2, 32.4]), while those with CC or TC genotypes had a median BDR of 10.9 (IQR = [5.0, 22.2]). SPATS2L mRNA knockdown resulted in increased ß(2)-adrenergic receptor levels. Our results suggest that SPATS2L may be an important regulator of ß(2)-adrenergic receptor down-regulation and that there is promise in gaining a better understanding of the biological mechanisms of differential response to ß(2)-agonists through GWAS.

Identification of a chronic obstructive pulmonary disease genetic determinant that regulates HHIP.

Multiple intergenic single-nucleotide polymorphisms (SNPs) near hedgehog interacting protein (HHIP) on chromosome 4q31 have been strongly associated with pulmonary function levels and moderate-to-severe chronic obstructive pulmonary disease (COPD). However, whether the effects of variants in this region are related to HHIP or another gene has not been proven. We confirmed genetic association of SNPs in the 4q31 COPD genome-wide association study (GWAS) region in a Polish cohort containing severe COPD cases and healthy smoking controls (P = 0.001 to 0.002). We found that HHIP expression at both mRNA and protein levels is reduced in COPD lung tissues. We identified a genomic region located ∼85 kb upstream of HHIP which contains a subset of associated SNPs, interacts with the HHIP promoter through a chromatin loop and functions as an HHIP enhancer. The COPD risk haplotype of two SNPs within this enhancer region (rs6537296A and rs1542725C) was associated with statistically significant reductions in HHIP promoter activity. Moreover, rs1542725 demonstrates differential binding to the transcription factor Sp3; the COPD-associated allele exhibits increased Sp3 binding, which is consistent with Sp3's usual function as a transcriptional repressor. Thus, increased Sp3 binding at a functional SNP within the chromosome 4q31 COPD GWAS locus leads to reduced HHIP expression and increased susceptibility to COPD through distal transcriptional regulation. Together, our findings reveal one mechanism through which SNPs upstream of the HHIP gene modulate the expression of HHIP and functionally implicate reduced HHIP gene expression in the pathogenesis of COPD.

Genome-wide association study reveals class I MHC-restricted T cell-associated molecule gene (CRTAM) variants interact with vitamin D levels to affect asthma exacerbations.

BACKGROUND: It has recently been shown that vitamin D deficiency can increase asthma development and severity and that variations in vitamin D receptor genes are associated with asthma susceptibility. OBJECTIVE: We sought to find genetic factors that might interact with vitamin D levels to affect the risk of asthma exacerbation. METHODS: We conducted a genome-wide study of gene-vitamin D interaction on asthma exacerbations using population-based and family-based approaches on 403 subjects and trios from the Childhood Asthma Management Program. Twenty-three polymorphisms with significant interactions were studied in a replication analysis in 584 children from a Costa Rican cohort. RESULTS: We identified 3 common variants in the class I MHC-restricted T cell-associated molecule gene (CRTAM) that were associated with an increased rate of asthma exacerbations based on the presence of a low circulating vitamin D level. These results were replicated in a second independent population (unadjusted combined interaction, P = .00028-.00097; combined odds ratio, 3.28-5.38). One variant, rs2272094, is a nonsynonymous coding polymorphism of CRTAM. Functional studies on cell lines confirmed the interaction of vitamin D and rs2272094 on CRTAM expression. CRTAM is highly expressed in activated human CD8(+) and natural killer T cells, both of which have been implicated in asthmatic patients. CONCLUSION: The findings highlight an important gene-environment interaction that elucidates the role of vitamin D and CD8(+) and natural killer T cells in asthma exacerbation in a genome-wide gene-environment interaction study that has been replicated in an independent population. The results suggest the potential importance of maintaining adequate vitamin D levels in subsets of high-risk asthmatic patients.

A genome-wide association study of COPD identifies a susceptibility locus on chromosome 19q13.

The genetic risk factors for chronic obstructive pulmonary disease (COPD) are still largely unknown. To date, genome-wide association studies (GWASs) of limited size have identified several novel risk loci for COPD at CHRNA3/CHRNA5/IREB2, HHIP and FAM13A; additional loci may be identified through larger studies. We performed a GWAS using a total of 3499 cases and 1922 control subjects from four cohorts: the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE); the Normative Aging Study (NAS) and National Emphysema Treatment Trial (NETT); Bergen, Norway (GenKOLS); and the COPDGene study. Genotyping was performed on Illumina platforms with additional markers imputed using 1000 Genomes data; results were summarized using fixed-effect meta-analysis. We identified a new genome-wide significant locus on chromosome 19q13 (rs7937, OR = 0.74, P = 2.9 × 10(-9)). Genotyping this single nucleotide polymorphism (SNP) and another nearby SNP in linkage disequilibrium (rs2604894) in 2859 subjects from the family-based International COPD Genetics Network study (ICGN) demonstrated supportive evidence for association for COPD (P = 0.28 and 0.11 for rs7937 and rs2604894), pre-bronchodilator FEV(1) (P = 0.08 and 0.04) and severe (GOLD 3&4) COPD (P = 0.09 and 0.017). This region includes RAB4B, EGLN2, MIA and CYP2A6, and has previously been identified in association with cigarette smoking behavior.

SOX5 is a candidate gene for chronic obstructive pulmonary disease susceptibility and is necessary for lung development.

RATIONALE: Chromosome 12p has been linked to chronic obstructive pulmonary disease (COPD) in the Boston Early-Onset COPD Study (BEOCOPD), but a susceptibility gene in that region has not been identified. OBJECTIVES: We used high-density single-nucleotide polymorphism (SNP) mapping to implicate a COPD susceptibility gene and an animal model to determine the potential role of SOX5 in lung development and COPD. METHODS: On chromosome 12p, we genotyped 1,387 SNPs in 386 COPD cases from the National Emphysema Treatment Trial and 424 control smokers from the Normative Aging Study. SNPs with significant associations were then tested in the BEOCOPD study and the International COPD Genetics Network. Based on the human results, we assessed histology and gene expression in the lungs of Sox5(-/-) mice. MEASUREMENTS AND MAIN RESULTS: In the case-control analysis, 27 SNPs were significant at P ≤ 0.01. The most significant SNP in the BEOCOPD replication was rs11046966 (National Emphysema Treatment Trial-Normative Aging Study P = 6.0 × 10(-4), BEOCOPD P = 1.5 × 10(-5), combined P = 1.7 × 10(-7)), located 3' to the gene SOX5. Association with rs11046966 was not replicated in the International COPD Genetics Network. Sox5(-/-) mice showed abnormal lung development, with a delay in maturation before the saccular stage, as early as E16.5. Lung pathology in Sox5(-/-) lungs was associated with a decrease in fibronectin expression, an extracellular matrix component critical for branching morphogenesis. CONCLUSIONS: Genetic variation in the transcription factor SOX5 is associated with COPD susceptibility. A mouse model suggests that the effect may be due, in part, to its effects on lung development and/or repair processes.

Genome-wide association analysis of body mass in chronic obstructive pulmonary disease.

Cachexia, whether assessed by body mass index (BMI) or fat-free mass index (FFMI), affects a significant proportion of patients with chronic obstructive pulmonary disease (COPD), and is an independent risk factor for increased mortality, increased emphysema, and more severe airflow obstruction. The variable development of cachexia among patients with COPD suggests a role for genetic susceptibility. The objective of the present study was to determine genetic susceptibility loci involved in the development of low BMI and FFMI in subjects with COPD. A genome-wide association study (GWAS) of BMI was conducted in three independent cohorts of European descent with Global Initiative for Chronic Obstructive Lung Disease stage II or higher COPD: Evaluation of COPD Longitudinally to Identify Predictive Surrogate End-Points (ECLIPSE; n = 1,734); Norway-Bergen cohort (n = 851); and a subset of subjects from the National Emphysema Treatment Trial (NETT; n = 365). A genome-wide association of FFMI was conducted in two of the cohorts (ECLIPSE and Norway). In the combined analyses, a significant association was found between rs8050136, located in the first intron of the fat mass and obesity-associated (FTO) gene, and BMI (P = 4.97 × 10(-7)) and FFMI (P = 1.19 × 10(-7)). We replicated the association in a fourth, independent cohort consisting of 502 subjects with COPD from COPDGene (P = 6 × 10(-3)). Within the largest contributing cohort of our analysis, lung function, as assessed by forced expiratory volume at 1 second, varied significantly by FTO genotype. Our analysis suggests a potential role for the FTO locus in the determination of anthropomorphic measures associated with COPD.

Mapping of numerous disease-associated expression polymorphisms in primary peripheral blood CD4+ lymphocytes.

Genome-wide association studies of human gene expression promise to identify functional regulatory genetic variation that contributes to phenotypic diversity. However, it is unclear how useful this approach will be for the identification of disease-susceptibility variants. We generated gene expression profiles for 22 184 mRNA transcripts using RNA derived from peripheral blood CD4+ lymphocytes, and genome-wide genotype data for 516 512 autosomal markers in 200 subjects. We screened for cis-acting variants by testing variants mapping within 50 kb of expressed transcripts for association with transcript abundance using generalized linear models. Significant associations were identified for 1585 genes at a false discovery rate of 0.05 (corresponding to P-values ranging from 1 × 10(-91) to 7 × 10(-4)). Importantly, we identified evidence of regulatory variation for 119 previously mapped disease genes, including 24 examples where the variant with the strongest evidence of disease-association demonstrates strong association with specific transcript abundance. The prevalence of cis-acting variants among disease-associated genes was 63% higher than the genome-wide rate in our data set (P = 6.41 × 10(-6)), and although many of the implicated loci were associated with immune-related diseases (including asthma, connective tissue disorders and inflammatory bowel disease), associations with genes implicated in non-immune-related diseases including lipid profiles, anthropomorphic measurements, cancer and neurologic disease were also observed. Genetic variants that confer inter-individual differences in gene expression represent an important subset of variants that contribute to disease susceptibility. Population-based integrative genetic approaches can help identify such variation and enhance our understanding of the genetic basis of complex traits.

Analyses of shared genetic factors between asthma and obesity in children.

BACKGROUND: Epidemiologic studies consistently show associations between asthma and obesity. Shared genetics might account for this association. OBJECTIVE: We sought to identify genetic variants associated with both asthma and obesity. METHODS: On the basis of a literature search, we identified genes from (1) genome-wide association studies (GWASs) of body mass index (BMI; n = 17 genes), (2) GWASs of asthma (n = 14), and (3) candidate gene studies of BMI and asthma (n = 7). We used GWAS data from the Childhood Asthma Management Program to analyze associations between single nucleotide polymorphisms (SNPs) in these genes and asthma (n = 359 subjects) and BMI (n = 537). RESULTS: One top BMI GWAS SNP from the literature, rs10938397 near glucosamine-6-phosphate deaminase 2 (GNPDA2), was associated with both BMI (P = 4 x 10(-4)) and asthma (P = .03). Of the top asthma GWAS SNPs and the candidate gene SNPs, none was found to be associated with both BMI and asthma. Gene-based analyses that included all available SNPs in each gene found associations (P < .05) with both phenotypes for several genes: neuronal growth regulator 1 (NEGR1); roundabout, axon guidance receptor, homolog 1 (ROBO1); diacylglycerol kinase, gamma (DGKG); Fas apoptotic inhibitory molecule 2 (FAIM2); fat mass and obesity associated (FTO); and carbohydrate (N-acetylgalactosamine 4-0) sulfotransferase 8 (CHST8) among the BMI GWAS genes; interleukin 1 receptor-like 1 / interleukin 18 receptor 1 (IL1RL1/IL18R1), dipeptidyl-peptidase 10 (DPP10), phosphodiesterase 4D (PDE4D), V-myb myeloblastosis viral oncogene homolog (MYB), PDE10A, IL33, and especially protein tyrosine phosphatase, receptor type D (PTPRD) among the asthma GWAS genes; and protein kinase C, alpha (PRKCA) among the BMI and asthma candidate genes. CONCLUSIONS: SNPs within several genes showed associations to BMI and asthma at a genetic level, but none of these associations were significant after correction for multiple testing. Our analysis of known candidate genes reveals some evidence for shared genetics between asthma and obesity, but other shared genetic determinants are likely to be identified in novel loci.

Fungal exposure modulates the effect of polymorphisms of chitinases on emergency department visits and hospitalizations.

RATIONALE: Chitinases are enzymes that cleave chitin, which is present in fungal cells. Two types of human chitinases, chitotriosidase and acidic mammalian chitinase, and the chitinase-like protein, YKL-40, seem to play an important role in asthma. We hypothesized that exposure to environmental fungi may modulate the effect of chitinases in individuals with asthma. OBJECTIVES: To explore whether interactions between high fungal exposure and common genetic variants in the two chitinases in humans, CHIT1 and CHIA, and the chitinase 3-like 1 gene, CHI3L1, are associated with severe asthma exacerbations and other asthma-related outcomes. METHODS: Forty-eight single nucleotide polymorphisms (SNPs) in CHIT1, CHIA, and CHI3L1 and one CHIT1 duplication were genotyped in 395 subjects and their parents as part of the Childhood Asthma Management Program. Household levels of mold (an index of fungal exposure) were determined on house dust samples. We conducted family-based association tests with gene-environment interactions. Our outcome was severe exacerbation, defined as emergency department visits and hospitalizations from asthma over a 4-year period, and our secondary outcomes included indices of lung function and allergy-related phenotypes. MEASUREMENTS AND MAIN RESULTS: Of the 395 subjects who had mold levels at randomization, 24% (95 subjects) had levels that were greater than 25,000 units per gram of house dust (high mold exposure). High mold exposure significantly modified the relation between three SNPs in CHIT1 (rs2486953, rs4950936, and rs1417149) and severe exacerbations (P for interaction 0.0010 for rs2486953, 0.0008 for rs4950936, and 0.0005 for rs1417149). High mold exposure did not significantly modify the relationship between any of the other variants and outcomes. CONCLUSIONS: Environmental exposure to fungi, modifies the effect of CHIT1 SNPs on severe asthma exacerbations.

Multistudy fine mapping of chromosome 2q identifies XRCC5 as a chronic obstructive pulmonary disease susceptibility gene.

RATIONALE: Several family-based studies have identified genetic linkage for lung function and airflow obstruction to chromosome 2q. OBJECTIVES: We hypothesized that merging results of high-resolution single nucleotide polymorphism (SNP) mapping in four separate populations would lead to the identification of chronic obstructive pulmonary disease (COPD) susceptibility genes on chromosome 2q. METHODS: Within the chromosome 2q linkage region, 2,843 SNPs were genotyped in 806 COPD cases and 779 control subjects from Norway, and 2,484 SNPs were genotyped in 309 patients with severe COPD from the National Emphysema Treatment Trial and 330 community control subjects. Significant associations from the combined results across the two case-control studies were followed up in 1,839 individuals from 603 families from the International COPD Genetics Network (ICGN) and in 949 individuals from 127 families in the Boston Early-Onset COPD Study. MEASUREMENTS AND MAIN RESULTS: Merging the results of the two case-control analyses, 14 of the 790 overlapping SNPs had a combined P < 0.01. Two of these 14 SNPs were consistently associated with COPD in the ICGN families. The association with one SNP, located in the gene XRCC5, was replicated in the Boston Early-Onset COPD Study, with a combined P = 2.51 x 10(-5) across the four studies, which remains significant when adjusted for multiple testing (P = 0.02). Genotype imputation confirmed the association with SNPs in XRCC5. CONCLUSIONS: By combining data from COPD genetic association studies conducted in four independent patient samples, we have identified XRCC5, an ATP-dependent DNA helicase, as a potential COPD susceptibility gene.

Variants in FAM13A are associated with chronic obstructive pulmonary disease.

We performed a genome-wide association study for chronic obstructive pulmonary disease (COPD) in three population cohorts, including 2,940 cases and 1,380 controls who were current or former smokers with normal lung function. We identified a new susceptibility locus at 4q22.1 in FAM13A and replicated this association in one case-control group (n = 1,006) and two family-based cohorts (n = 3,808) (rs7671167, combined P = 1.2 x 10(-11), combined odds ratio in case-control studies 0.76, 95% confidence interval 0.69-0.83).

MMP12, lung function, and COPD in high-risk populations.

BACKGROUND: Genetic variants influencing lung function in children and adults may ultimately lead to the development of chronic obstructive pulmonary disease (COPD), particularly in high-risk groups. METHODS: We tested for an association between single-nucleotide polymorphisms (SNPs) in the gene encoding matrix metalloproteinase 12 (MMP12) and a measure of lung function (prebronchodilator forced expiratory volume in 1 second [FEV(1)]) in more than 8300 subjects in seven cohorts that included children and adults. Within the Normative Aging Study (NAS), a cohort of initially healthy adult men, we tested for an association between SNPs that were associated with FEV(1) and the time to the onset of COPD. We then examined the relationship between MMP12 SNPs and COPD in two cohorts of adults with COPD or at risk for COPD. RESULTS: The minor allele (G) of a functional variant in the promoter region of MMP12 (rs2276109 [-82A-->G]) was positively associated with FEV(1) in a combined analysis of children with asthma and adult former and current smokers in all cohorts (P=2x10(-6)). This allele was also associated with a reduced risk of the onset of COPD in the NAS cohort (hazard ratio, 0.65; 95% confidence interval [CI], 0.46 to 0.92; P=0.02) and with a reduced risk of COPD in a cohort of smokers (odds ratio, 0.63; 95% CI, 0.45 to 0.88; P=0.005) and among participants in a family-based study of early-onset COPD (P=0.006). CONCLUSIONS: The minor allele of a SNP in MMP12 (rs2276109) is associated with a positive effect on lung function in children with asthma and in adults who smoke. This allele is also associated with a reduced risk of COPD in adult smokers.

PRKCA: a positional candidate gene for body mass index and asthma.

Asthma incidence and prevalence are higher in obese individuals. A potential mechanistic basis for this relationship is pleiotropy. We hypothesized that significant linkage and candidate-gene association would be found for body mass index (BMI) in a population ascertained on asthma affection status. Linkage analysis for BMI was performed on 657 subjects in eight Costa Rican families enrolled in a study of asthma. Family-based association studies were conducted for BMI with SNPs within a positional candidate gene, PRKCA. SNPs within PRKCA were also tested for association with asthma. Association studies were conducted in 415 Costa Rican parent-child trios and 493 trios participating in the Childhood Asthma Management Program (CAMP). Although only modest evidence of linkage for BMI was obtained for the whole cohort, significant linkage was noted for BMI in females on chromosome 17q (peak LOD = 3.39). Four SNPs in a candidate gene in this region (PRKCA) had unadjusted association p values < 0.05 for BMI in both cohorts, with the joint p value for two SNPs remaining significant after adjustment for multiple comparisons (rs228883 and rs1005651, joint p values = 9.5 x 10(-)(5) and 5.6 x 10(-)(5)). Similarly, eight SNPs had unadjusted association p values < 0.05 for asthma in both populations, with one SNP remaining significant after adjustment for multiple comparisons (rs11079657, joint p value = 2.6 x 10(-)(5)). PRKCA is a pleiotropic locus that is associated with both BMI and asthma and that has been identified via linkage analysis of BMI in a population ascertained on asthma.

Platform dependence of inference on gene-wise and gene-set involvement in human lung development.

BACKGROUND: With the recent development of microarray technologies, the comparability of gene expression data obtained from different platforms poses an important problem. We evaluated two widely used platforms, Affymetrix U133 Plus 2.0 and the Illumina HumanRef-8 v2 Expression Bead Chips, for comparability in a biological system in which changes may be subtle, namely fetal lung tissue as a function of gestational age. RESULTS: We performed the comparison via sequence-based probe matching between the two platforms. "Significance grouping" was defined as a measure of comparability. Using both expression correlation and significance grouping as measures of comparability, we demonstrated that despite overall cross-platform differences at the single gene level, increased correlation between the two platforms was found in genes with higher expression level, higher probe overlap, and lower p-value. We also demonstrated that biological function as determined via KEGG pathways or GO categories is more consistent across platforms than single gene analysis. CONCLUSION: We conclude that while the comparability of the platforms at the single gene level may be increased by increasing sample size, they are highly comparable ontologically even for subtle differences in a relatively small sample size. Biologically relevant inference should therefore be reproducible across laboratories using different platforms.

Assessing the reproducibility of asthma candidate gene associations, using genome-wide data.

RATIONALE: Association studies have implicated many genes in asthma pathogenesis, with replicated associations between single-nucleotide polymorphisms (SNPs) and asthma reported for more than 30 genes. Genome-wide genotyping enables simultaneous evaluation of most of this variation, and facilitates more comprehensive analysis of other common genetic variation around these candidate genes for association with asthma. OBJECTIVES: To use available genome-wide genotypic data to assess the reproducibility of previously reported associations with asthma and to evaluate the contribution of additional common genetic variation surrounding these loci to asthma susceptibility. METHODS: Illumina Human Hap 550Kv3 BeadChip (Illumina, San Diego, CA) SNP arrays were genotyped in 422 nuclear families participating in the Childhood Asthma Management Program. Genes with at least one SNP demonstrating prior association with asthma in two or more populations were tested for evidence of association with asthma, using family-based association testing. MEASUREMENTS AND MAIN RESULTS: We identified 39 candidate genes from the literature, using prespecified criteria. Of the 160 SNPs previously genotyped in these 39 genes, 10 SNPs in 6 genes were significantly associated with asthma (including the first independent replication for asthma-associated integrin beta(3) [ITGB3]). Evaluation of 619 additional common variants included in the Illumina 550K array revealed additional evidence of asthma association for 15 genes, although none were significant after adjustment for multiple comparisons. CONCLUSIONS: We replicated asthma associations for a minority of candidate genes. Pooling genome-wide association study results from multiple studies will increase the power to appreciate marginal effects of genes and further clarify which candidates are true "asthma genes."

Transforming growth factor-beta receptor-3 is associated with pulmonary emphysema.

Chronic obstructive pulmonary disease (COPD) is a heterogeneous syndrome, including emphysema and airway disease. Phenotypes defined on the basis of chest computed tomography (CT) may decrease disease heterogeneity and aid in the identification of candidate genes for COPD subtypes. To identify these genes, we performed genome-wide linkage analysis in extended pedigrees from the Boston Early-Onset COPD Study, stratified by emphysema status (defined by chest CT scans) of the probands, followed by genetic association analysis of positional candidate genes. A region on chromosome 1p showed strong evidence of linkage to lung function traits in families of emphysema-predominant probands in the stratified analysis (LOD score = 2.99 in families of emphysema-predominant probands versus 1.98 in all families). Association analysis in 949 individuals from 127 early-onset COPD pedigrees revealed association for COPD-related traits with an intronic single-nucleotide polymorphism (SNP) in transforming growth factor-beta receptor-3 (TGFBR3) (P = 0.005). This SNP was significantly associated with COPD affection status comparing 389 cases from the National Emphysema Treatment Trial to 472 control smokers (P = 0.04), and with FEV(1) (P = 0.004) and CT emphysema (P = 0.05) in 3,117 subjects from the International COPD Genetics Network. Gene-level replication of association with lung function was seen in 427 patients with COPD from the Lung Health Study. In conclusion, stratified linkage analysis followed by association testing identified TGFBR3 (betaglycan) as a potential susceptibility gene for COPD. Published human microarray and murine linkage studies have also demonstrated the importance of TGFBR3 in emphysema and lung function, and our group and others have previously found association of COPD-related traits with TGFB1, a ligand for TGFBR3.

Analysis of exonic elastin variants in severe, early-onset chronic obstructive pulmonary disease.

The destruction of elastic fibers has been implicated in the pathogenesis of chronic obstructive pulmonary disease (COPD). Emphysema has been described in autosomal dominant cutis laxa, which can be caused by mutations in the elastin gene. Previously, a rare functional mutation in the terminal exon of elastin was found in a case of severe, early-onset COPD. To test the hypothesis that other similar elastin mutations may predispose to COPD, we screened 90 probands from the Boston Early-Onset COPD Study and 90 smoking control subjects from the Normative Aging Study for mutations in elastin exons using high-resolution DNA melt analysis followed by resequencing. Rare nonsynonymous single-nucleotide polymorphisms (SNPs) seen only in cases were examined for segregation with airflow obstruction within pedigrees. Common nonsynonymous SNPs were tested for association with COPD in a family-based analysis of 949 subjects from the Boston Early-Onset COPD Study, and in a case-control analysis in 389 COPD cases from the National Emphysema Treatment Trial and 472 control subjects from the Normative Aging Study. Of 28 elastin variants found, 3 were nonsynonymous SNPs found only in cases. The previously described Gly773Asp mutation was found in another proband. The other two SNPs did not clearly segregate with COPD within families. Two common nonsynonymous SNPs did not demonstrate significant associations in either a family-based or case-control analysis. Exonic SNPs in the elastin gene do not appear to be common risk factors for severe COPD.